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Abstract 

While built on a sound physical foundation, isentropic, one-dimensional models generally used to analyze the dynamics of dilute two-phase 
(feed-powder particles suspended in a carrier gas) flow during the cold-gas dynamic-spray process, require the use of numerical procedures 
to obtain solutions for the governing equations. Numerical solutions, unfortunately, do not enable an easy establishment of the relationships 
between the gas, process and feed-powder parameters on one side and the gas and the particle velocities at the nozzle exit and the particle 
impact velocity, on the other. Analytical solutions for the governing equations in the limits of small and large relative particle/gas velocities 
and a multiple non-linear regression analysis are used, in the present work, to develop analytical functions which can be used to compute the 
gas and the particle exit velocities and the particle impact velocity for a given set of the gas, process, and feed-powder parameters. The results 
obtained using the analytical functions are found to be in a very good agreement with their numerical and experimental counterparts. 

© 2003 Elsevier B.V. All rights reserved. 
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1. Introduction 

The cold-gas dynamic-spray process, often referred to 
as simply “cold spray,” is a high-rate material deposition 
process in which fine, solid powder particles (generally 
1-50 ptm in diameter) are accelerated in a supersonic jet of 
compressed (carrier) gas to velocities in a range between 
500 and 1000 m/s. As the solid particles impact the target 
surface, they undergo plastic deformation and bond to the 
surface, rapidly building up a layer of deposited material. 
Cold spray as a coating technology was initially developed 
in the mid-1980s at the Institute for Theoretical and Ap¬ 
plied Mechanics of the Siberian Division of the Russian 
Academy of Science in Novosibirsk [1,2]. The Russian sci¬ 
entists successfully deposited a wide range of pure metals, 
metallic alloys, polymers, and composites onto a variety of 
substrate materials. In addition, they demonstrated that very 
high coating deposition rates on the order of 5 m 2 /min (~ 
300ft 2 /min) are attainable using the cold-spray process. 
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A simple schematic of a typical cold-spray device is 
shown in Fig. 1 . Compressed gas of an inlet pressure on the 
order of 30 bar (500 psi) enters the device and flows through 
a converging/diverging nozzle to attain a supersonic ve¬ 
locity. The solid powder particles are metered into the gas 
flow upstream of the converging section of the nozzle and 
are accelerated by the rapidly expanding gas. To achieve 
higher gas flow velocities in the nozzle, the compressed 
gas is often preheated. However, while preheat tempera¬ 
tures as high as 900 K are sometimes used, due to the fact 
that the contact time of spray particles with the hot gas is 
quite short and that the gas rapidly cools as it expands in 
the diverging section of the nozzle, the temperature of the 
particles remains substantially below the initial gas preheat 
temperature and, hence, below the melting temperature of 
the powder material. 

The actual mechanism by which the solid particles deform 
and bond during cold spray is still not well understood. The 
prevailing theory for cold-spray bonding postulates that, dur¬ 
ing impact, the solid particles undergo plastic deformation, 
disrupt thin (oxide) surface films and, in turn, achieve inti¬ 
mate conformal contact with the target surface. The intimate 
conformal contact combined with high contact pressures 
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Nomenclature 

A 

nozzle cross-section area 

C 

specific heat 

Cd 

drag coefficient 

D 

drag force 

h 

convective heat transfer coefficient 

L 

thickness of the stagnant region 

m 

mass 

m 

mass flow rate 

p 

pressure 

R 

gas constant 

Re 

Reynolds number 

S 

molecular speed ratio 

T 

temperature 

V 

velocity 

X 

axial distance from the nozzle throat 

p 

viscosity 

Y 

specific heat ratio 

P 

density 

Subscripts 

0 

stagnation quantity 

e 

nozzle-exit related quantity 

p 

particle quantity 

s 

shock related quantity 

Superscripts 

* 

nozzle throat quantity 

impact 

impact related quantity 

opt 

optimal quantity 

St 

stagnant region related quantity 


promotes bonding. This theory is supported by a number 
of experimental findings such as: (a) a wide range of duc¬ 
tile (metallic and polymeric) materials can be successfully 
cold-sprayed while non-ductile materials such as ceramics 
can be deposited only if they are co-cold-sprayed with a 
ductile (matrix) material [2]; (b) the mean deposition parti¬ 
cle velocity should exceed a minimum (material-dependent) 
critical velocity to achieve deposition which suggests that 


Powder Feeder 



Substrate 


Fig. 1. Schematic of a typical cold-gas dynamic-spray system. 


sufficient kinetic energy must be available to plastically 
deform the solid material and/or disrupt the surface film 
[3]; and (c) the particle kinetic energy at impact is typically 
significantly lower than the energy required to melt the par¬ 
ticle suggesting that the deposition mechanism is primarily, 
or perhaps entirely, a solid-state process [4-8]. 

As the cold-spray process does not normally involve the 
use of a high-temperature heat source, it generally offers 
a number of advantages over the thermal-spray material 
deposition technologies such as high velocity oxy-fuel, 
detonation gun, plasma spray, and arc spray. Among these 
advantages, the most important appear to be: (a) the amount 
of heat delivered to the coated part is relatively small so 
that microstructural changes in the substrate material are 
minimal or nonexistent; (b) due to the absence of in-flight 
oxidation and other chemical reactions, thermally- and 
oxygen-sensitive depositing materials (e.g. copper or ti¬ 
tanium) can be cold sprayed without significant material 
degradation; (c) nanophase, intermetallic and amorphous 
materials, which are not amenable to conventional ther¬ 
mal spray processes (due to a major degradation of the 
depositing material), can be cold sprayed; (d) formation 
of the embrittling phases is generally avoided; (e) macro- 
and micro-segregations of the alloying elements during so¬ 
lidification which accompany conventional thermal spray 
techniques and can considerably compromise materials 
properties do not occur during cold spraying. Consequently, 
attractive properties are retained in cold-sprayed bulk ma¬ 
terials; (f) “peening” effect of the impinging solid particles 
can give rise to potentially beneficial compressive residual 
stresses in cold-spray deposited materials [3] in contrast 
to the highly detrimental tensile residual stresses induced 
by solidification shrinkage accompanying the conventional 
thermal-spray processes; and (g) cold spray of the materials 
like copper, solder and polymeric coatings offers exciting 
new possibilities for cost-effective and environmentally- 
friendly alternatives to technologies such as electroplating, 
soldering and painting [9]. 

The solid-state particle-bonding mechanism discussed 
above suggests and experimental observations [9] con¬ 
firm that it is desirable to maximize the velocity at which 
feed-powder particles impact the target surface. While, in 
general, this can be accomplished by increasing the inlet 
pressure of the carried gas, for practical and economic rea¬ 
sons, it is desirable to maximize the particle impact velocity 
at a given level of the carrier-gas inlet pressure by properly 
selecting the type of the carrier gas, its inlet temperature 
and by optimizing the shape of the converging/diverging 
cold-spray nozzle. A detailed analysis of the effects that the 
type of the carrier gas, the inlet gas temperature and the 
shape of the cold-spray nozzle have on the impact velocity 
of the feed-powder particles is carried out by Dykhuizen and 
Smith [10] using an isentropic, one-dimensional gas-flow 
model. In our recent work [11], the analysis of Dykhuizen 
and Smith [10] was extended in order to include the effects 
of finite values of the particle velocity and the effect of 
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variability of the gas/particle drag coefficient. While the 
one-dimensional model of Dykhuizen and Smith [10] has 
been found to be quite successful, its numerical nature does 
not enable an easy establishment of the relationships be¬ 
tween the gas, process and feed-powder parameters on one 
side and the gas and the particle velocities at the nozzle exit 
and the velocity at which particles impact the substrate sur¬ 
face, on the other. The objective of the present paper is to 
use the results of the one-dimensional model by Dykhuizen 
and Smith [10] and a multiple regression analysis to de¬ 
velop analytical functions which can be used to compute 
the gas and the particle exit velocities for a given set of the 
gas, process and feed-powder parameters. 

The organization of the paper is as follows: A brief 
overview of the one-dimensional isoentropic gas flow model 
developed by Dykhuizen and Smith [10] utilized in the 
present work is presented in Section 2. The one-dimensional 
particle dynamics model also proposed by Dykhuizen and 
Smith [10] is briefly discussed in Section 3. The derivation 
of the analytical functions for gas and powder-particle ve¬ 
locities at the nozzle exit and for the impact particle velocity 
are presented in Sections 4.1-4.3, respectively. The key 
conclusions resulted from the present study are summarized 
in Section 5. 


2. Isentropic gas flow model 


In this section, a brief overview is given of the cold-spray 
gas-flow model developed by Dykhuizen and Smith [10]. 
The model considers a typical geometry of the cold-spray 
converging/diverging nozzle (Fig. 1) and involves a number 
of assumptions and simplifications such as: (a) the gas flow 
is assumed to be one dimensional and isentropic (adiabatic 
and frictionless); (b) the gas is treated as a perfect (ideal) 
gas; and (c) the constant-pressure and the constant-volume 
specific heats of the gas are assumed to be constant. 

The carrier gas flow is assumed to originate from a large 
chamber or duct where its velocity is zero and the pressure 
(referred to as the “stagnation” pressure) is P a and the tem¬ 
perature (referred to as the “total” gas temperature) is T 0 . The 
cold-spray process is furthermore assumed to be controlled 
by the user who can set the total temperature and the mass 
flow rate of the gas. The corresponding stagnation pressure 
can then be calculated using the following procedure: 

Using basic dynamic and thermodynamic relations for the 
compressible fluid flow, the gas temperature T* at the small¬ 
est cross-sectional area of the converging/diverging nozzle 
(referred to as the nozzle throat in the following) where the 
Mach number (the velocity of the gas divided by the local 
speed of sound) is unity, can be derived as [12]: 



(1-l-y-1/2) 


( 1 ) 


where y is the ratio of the constant-pressure and the 
constant-volume specific heats which is typically set to 1.66 


for monoatomic gases like helium and 1.4 for diatomic 
gases like nitrogen and oxygen, respectively. Since air is pri¬ 
marily a mixture of nitrogen and oxygen, it is considered as 
a diatomic gas. y takes on values smaller than 1.4 in multi¬ 
atom molecular gases but such gases are rarely used in the 
cold-spray process. 

The gas velocity at the nozzle throat is equal to the speed 
of sound and is defined as: 

( 2 ) 


where R is the gas constant (the universal gas constant di¬ 
vided by the gas molecular weight). It should be noted that 
the superscript * is used throughout this paper to denote the 
quantities at the nozzle throat, that is the gas quantities un¬ 
der the sonic conditions. Eq. (2) can be used to explain the 
experimental finding that the low molecular weight (and, 
hence, large R ), monoatomic (and, hence, large y and, in 
turn, high T *) helium is a better carrier gas than the high 
molecular weight, diatomic air since for the same total gas 
temperature T 0 , it is associated with a higher speed of sound. 

From the known (user selected) mass flow rate, m, the 
sonic gas density can be computed as: 


where A* is the (known) cross-sectional area of the nozzle 
throat. 

Next, using the ideal gas law, the gas pressure at the nozzle 
throat can be determined as: 


P* = p*RT* 


(4) 


Once the throat pressure P* is computed, the stagnation 
pressure P 0 can be calculated using the following isentropic 
relation: 



yKy- l) 


(5) 


After all the gas-dynamics quantities ( T *, v*, p*, and P*) are 
calculated at the nozzle throat, one can proceed to determine 
these quantities along the diverging section of the nozzle. 
Toward that end, the variation of one of the these quantities 
or the variation of the Mach number or the variation of the 
nozzle cross-sectional area along the diverging section of the 
nozzle must be specified. Dykhuizen and Smith [10] consid¬ 
ered the case when the variation of the nozzle cross-sectional 
area A is specified, which then allows determination of 
the corresponding Mach number M from the following 
equation: 


A 

A* 




Y/(2 y -l) 


( 6 ) 


Once the Mach number is determined at a given 
cross-sectional area of the diverging section of the nozzle, 
the remaining corresponding gas quantities (P, T, v, and p) 
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can be calculated using the following isentropic relation- 


ships: 


/>/.*( J ' +1 V' 1 ’”" 

(7) 

\2+(y-l)M 2 J 

T 0 


(l + (y-l/2)M 2 ) 

(8) 

u = My/yRT 

(9) 

Po 

P ~ (1 + (y- l)/ 2 )M 2 )V(y-i) 

(10) 


Eqs. (7)-(10) can be used to compute P, T, u, and p at 
the nozzle exit, if the given exit cross-sectional area A e is 
substituted for A in Eq. (6). However, it must be noted that P, 
T , V, and p defined in this way will reflect the true conditions 
of the gas at the nozzle exit only if a normal shock does 
not take place inside the nozzle. To determine if the normal 
shock will take place inside the nozzle, the ambient pressure 
should be compared with the following “shock” pressure: 


P s = P e 




( 11 ) 


where subscript e is used to denote the gas quantities at the 
nozzle exit. 

If the shock pressure P s is lower than the ambient pres¬ 
sure, a shock will occur inside the nozzle and the subsequent 
gas flow is subsonic so that the exit pressure is not given by 
Eq. (7), but it is rather equal to the ambient pressure. Un¬ 
der normal cold-spray operating conditions, the shock pres¬ 
sure is maintained above the ambient pressure so that no 
shock occurs inside the nozzle and the exit pressure is de¬ 
fined by Eq. (7). At the same time, the exit gas pressure, 
P e , is generally lower than the ambient pressure in effort 
to maximize the exit velocity of the gas (and thus, the av¬ 
erage velocity of the feed-powder particles, referred to as 
the particle velocity in the following). Under such condi¬ 
tions, the one-dimensional gas dynamics model developed 
by Dykhuizen and Smith [10] and briefly reviewed above 
can be used to analyze the gas flow inside the nozzle. As 
the gas leaves the nozzle, it slows down as the gas pressure 
tries to adjust to the ambient pressure. However, due to a rel¬ 
atively short nozzle-exit/substrate-surface standoff distance 
encountered in the cold-spray process, this decrease in the 
gas velocity is not expected to be significant. However, the 
impingement of the gas jet upon the substrate surface gives 
rise to the formation of a bow shock (stagnant) region in 
front of the substrate. The component of the gas velocity 
normal to the substrate surface in this region is quite small 
relative to the nozzle-exit gas and particle velocities. In ad¬ 
dition, the gas flow in this region is dominated by its lat¬ 
eral component. A more detailed discussion of the stagnant 
region and its effect of the particle impact velocity is pre¬ 
sented in Section 4.3. 


3. Particle dynamics model 


Dykhuizen and Smith [5] also analyzed the interactions 
of the carrier gas with the spray particles under the approx¬ 
imation of a dilute two-phase (gas + non-interacting solid 
particles) flow. The particle velocity, up, can, in this case, be 
determined by solving the following differential equation: 

C D A p p(n - Up) 2 


du n 


dt>n 


dr 


= m nil, 


p " p dx ~ 


( 12 ) 


where m p and A p are the average mass and the cross-sectional 
area of the particles, respectively, Cd the drag coefficient, 
t the time and x the axial distance traveled by the particle 
(measured from the nozzle throat). Under the condition 
of constant gas velocity, gas density and drag coefficient, 
Eq. (12) can be integrated to yield: 



v 

H- 

V ~ Up 




CdA p px 
2m p 


(13) 


It should be noted that Eq. (13) may generally not be valid 
since both the gas density and the drag coefficient may vary 
over the length of the nozzle [10,11]. 

A simple analysis of Eq. (12) shows that the ultimate 
particle velocity is equal to the gas velocity. Furthermore, 
examination of Eqs. (6), (8) and (9) indicates that the gas 
velocity within the nozzle depends on the total gas temper¬ 
ature and the nozzle geometry (i.e. the cross-sectional area 
at a given axial distance x), but not on the gas pressure (un¬ 
der the condition of a constant drag coefficient). However, 
Eqs. (9), (10) and (12) indicate that the initial particle ac¬ 
celeration (dup/d t at Up = 0) is linearly dependent on the 
stagnation pressure but independent of the total temperature. 
Thus, while the stagnation pressure does not affect the max¬ 
imum particle velocity, it has to be sufficiently high to en¬ 
sure that the spray-particles velocity will approach the gas 
velocity over a relatively short length of the diverging sec¬ 
tion of the nozzle. 


4. Results and discussion 

4.1. Gas velocity at the nozzle exit 

As discussed in the Section 2, the gas velocity at the 
nozzle exit can be computed using the one-dimensional nu¬ 
merical model of Dykhuizen and Smith [10]. However, this 
model entails the solution of a non-linear algebraic equa¬ 
tion, Eq. (6), and does not enable an easy access to the rela¬ 
tionships between the carrier-gas and the cold-spray process 
parameters on one side and the gas velocity at the nozzle 
exit, on the other. 

An attempt was made in this section to construct an ana¬ 
lytical function which can be used to compute the gas exit 
velocity for a given set of gas and process parameters. To¬ 
wards that end, Eq. (6) is solved for the Mach number M 
for different values of the nozzle expansion ratio, A/A*, and 
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Fig. 2. Variation of the gas Mach number with the nozzle expansion ratio 
and the gas specific-heat ration, y. 


for several values of the specific heat ratio, y. The results 
obtained are displayed as solid lines in the Fig. 2. Next, a 
non-linear least squares procedure is used to fit the M versus 
A!A* data for different y values to a function: 


M 



■ 


+ (l-*i) 


(14) 


where k\ and ki are function of the specific heat ratio, y. 
A non-linear polynomial regression analysis is next used to 
establish the y-dependence of k\ and k 2 as: 

ki = 218.0629- 243.5764k + 71.7925k 2 (15) 

k 2 = -0.122450 + 0.281300k (16) 


A comparison of the M versus At A * relation defined by 
Eq. (14), and the one predicted by the model of Dykhuizen 
and Smith [10], Eq. (6), is shown in Fig. 2. The M versus 
AlA* data based on Eq. (14) are shown as solid circles in 
Fig. 2. A reasonably good agreement is seen between the two 
sets of results particularly for high y values encountered in 
mono-atomic (e.g. helium) and di-atomic (e.g. air, nitrogen) 
gases commonly used in the cold-spray process. 

Next Eq. (8) is substituted in Eq. (9) to obtain: 

v — M I - yRT ° (17) 

V 1 + (K- 1/2)M 2 

Thus Eq. (17), in conjunction with Eqs. (14)—(16), can 
be used to compute the exit gas velocity v e as a function of 
the nozzle exit expansion ratio, (AJA*), and the carrier-gas 
(stagnation) properties: y, R, and 7 0 . 


4.2. Particle velocity at the nozzle exit 


The non-linear first-order ordinary differential equation, 
Eq. (12), does not have an analytical (closed form) solu¬ 


tion over the entire range of nozzle lengths, 0 < x < 00 . 
Therefore, a numerical integration routine must be used to 
solve Eq. (12). While an integration procedure can be read¬ 
ily implemented, the numerical solution for the u p versus x 
relationship obtained is not very attractive, since it does not 
offer an insight into the effect of various carrier gas, pro¬ 
cess, and feed-powder parameters. To overcome this limita¬ 
tion, analytical solutions for Eq. (12) are first sought in two 
distinct limits: (a) i> p <5jC i> and (b) u p v. 

When the particle velocity is small in comparison with the 
gas velocity, and the changes in gas density are negligible 
(p — constant), Eq. (12) can be solved to yield: 


v 0 = 


1 3 pCnx 

PpDp 


(18) 


where the gas density is given in terms of the stagnation gas 
density, p 0 , and the Mach number, M. via Eq. (10), while 
p 0 is related to the stagnation pressure, P 0 , and the total gas 
temperature, T 0 , via Eq. (4). As will be demonstrated later 
in this section, by comparing model predictions with their 
experimental counterparts the p = constant assumption used 
in the derivation of Eq. (18) is found to be justified. 

When the particle velocity is comparable with the gas 
velocity, (i.e. when the relative particle/gas Reynolds number 
is small), the particle velocity gradient is defined as: 

dl>n 18u(l>— 1> D ) 

df=- 5 *r- (19) 


where p is the gas viscosity which enters Eq. (19) through 
a relationship between the drag coefficient and the relative 
Reynolds number. Integration of Eq. (19) yields: 


v p = v 


Wo 



/ 18 R-x _ \ 
\ D pPp v / 



( 20 ) 


where Wo is the principal branch of the Lambert W function. 

The Lambert W function is defined as the function that 
solves the equation W(z)& w ^ = z where z is a real or a 
complex number. Since this equation always has an infi¬ 
nite number of solutions, most of them complex, W(z ) is a 
multi-valued function. The different possible solutions are 
labeled by an integer variable k called the branch of W(z). 
Thus, the proper way to talk about the solutions of the equa¬ 
tion W(z)e^ J l = Z is to say that they are: Wt(z) where k = 
0, ±1, ±2, etc. When z is a real number and z < — 1/e, the 
solution to the equation W(z)e w ^ = zis multi-valued and 
complex. When — 1 /e < z < 0 there are two possible real 
values of W(z). The branch satisfying W(z) > — 1 is denoted 
as Wo(z) and is called the principal branch. The branch sat¬ 
isfying W(z) < — 1 is denoted as W_i(z). When z is real 
and z > 0, the solution to the equation W(z)e w ^ = z is 
single-valued and real and belongs to the principal branch, 
Wo(z). The real solution of the equation W(z.)c' v(z] = z are 
displayed in Fig. 3. It is the portion of the principal branch 
for — 1 /e < z < 0 which is of interest in the present work 
since for 0 < x < 00 the magnitude of the Lambert function 
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Fig. 3. Real branches of the Lambert W function. 


should change between — 1 and 0 (to ensure that v p changes 
between 0 and v) and for x= 0 , the argument of the function 
W{)(z) in Eq. (20) is equal to —1/e while for x becoming 
infinite the argument of the function Wo(z) in Eq. (20) is 
zero. The Lambert W function is implemented in many com¬ 
mercial general-purpose computing systems such as Matlab 
(used in the present work), Maple, Macsyma, and Mathe- 
matica. For — 1 /e < z < 0, of interest in the present work, 
Wq(z) can be approximated using the following nonlinear 
function: 


-l -9 o 

m i = — e »?2 = —e m 3 = — Le 

.si = \/2 S 2 = 2 si — 3 S 3 = 4 — 3j2 S 4 = si — 2 


Wq(z) = a/ 



«/ 



a(s 2 P + s 3 ) \ 
s 4 (a + P) ) 


- 1 ( 21 ) 


Particle velocities normalized by the gas velocity, v p /v, 
as a function of the length of the diverging section of the 
nozzle, x, as predicted by Eqs. (18) and (20), are plotted in 
Fig. 4 (curves labeled “Eqs. (18) and (20)”). Eq. (18) is used 
at low values of x where the v p /v values are small while 
Eq. (20) is used at large values of x where the v p /v values 
are near unity. The following (typical) cold-spray process¬ 
ing conditions are used to generate the results presented in 
Fig. 4: carrier gas—helium, the total gas temperature T a = 
600 K, the stagnation pressure P 0 = 22 bar, gas viscosity 
jji = 1.9 x 10~ ? Ns/m 2 ; the mean particle diameter Dp = 
10 pun and the particle material density p p = 2.7 g/cm 3 . The 
variation of v p /v with x obtained using the one-dimensional 
numerical model of Dykhuizen and Smith [10] is also shown 
in Fig. 4 (curve labeled “Numerical”). It is seen that the an¬ 
alytical solution given by Eq. (18) agrees quite well with the 
numerical solution in the regions of its validity. The same 
holds for the analytical solution given by Eq. (20), but at 
very large x values. Since such x values are not of practical 


interest in the cold-gas dynamic-spray process, they are not 
shown in Fig. 4. 

To construct an analytical v p /v versus x function for the 
entire 0 < v p /v < lrange, it is next assumed that the v p /v 
versus x relationship for the intermediate v p /v values can 
be obtained as a linear combination of the two relations 
predicted by Eqs. (18) and (20). However, before such linear 
combination is constructed, the two v p /v versus x functional 
relations predicted by Eqs. (18) and (20) are respectively 
transformed into the following exponential forms: 

( l ^) = _e-V 3 P c WPpD P + , <<); ( 22 ) 

V V / 18 

W = + 1 v ^v (23) 

V v /20 


Eqs. (22) and (23) are next combined as: 


_T 

v 


= «>(-) +(!-«)(-) 

V V / 18 \|j / 20 


(24) 


where w and 1 — w are the weighting factors. In princi¬ 
ple, optimal values of the weighting factors can be obtained 
by minimizing deviations of the v p /v versus x relation pre¬ 
dicted by Eq. (24) from the numerical solution. This is not 
done in the present work to avoid unnecessary complexities 
associated with the dependence of such optimal values on 
various gas, process and feed-powder parameters. A value 
of w = 0.5 is found to give a reasonably good agreement 
between the u p /v versus x relation predicted by Eq. (24) and 
its numerical counterpart for a wide range of gas, process 
and feed-powder parameters. A typical level of agreement 
between the v p /v versus x relation predicted by Eq. (24) (de¬ 
noted as “Eq. (24)”) and its numerical counterpart (denoted 
as “Numerical”) is shown in Fig. 4. 



Fig. 4. Variation of the ratio of the particle velocity and the gas velocity 
at the nozzle exit with length of the nozzle. Please see text for details. 
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Thus, Eq. (24), in conjunction with Eqs. (22) and (23), 
can be used to compute the particle velocity at the nozzle 
exit for a given set of the gas parameters: y. /x, R, T 0 , and 
P 0 ; the process parameters: AM* and x\ and the feed-powder 
parameters: p p and D p . 

4.3. Particle impact velocity 

As discussed in Section 2, calculations of the velocity of 
the velocity at which feed-powder particles impact the sub¬ 
strate surface must include the effects associated with the 
impingement of the supersonic gas jet upon the substrate. Su¬ 
personic impinging jets have been the focus of research for 
over three decades because of their application in a number 
of diverse areas ranging from hovering of the short takeoff 
and vertical landing (STOVL) aircrafts to material process¬ 
ing (e.g. the cold-gas dynamic-spray process). While fluid 
dynamics and acoustic properties of the supersonic imping¬ 
ing jets vary considerably with the geometrical parameters 
of the nozzle, the distance between the nozzle exit and the 
substrate, the exit Mach number, etc., some common global 
features of the supersonic impinging jets can be identified. A 
schematic of the supersonic impinging jet flowfield is given 
in Fig. 5. It is seen that the impinging jet flowfield can be 
divided into three main regions [13]. The first region repre¬ 
sents the main jet column in which the flow is primarily in- 
viscid and contains expansion and compression shock waves 
for non-ideally expanding jets. The second region, generally 
referred as the impinging zone, involves the region of jet 
impingement onto the substrate. The impingement zone is 
characterized by large gradients which cause major changes 
in the local flow properties. A stagnation bubble containing 
re-circulating fluid with relatively low velocity is also de¬ 
picted within the impingement region in Fig. 5. The origin 
of this bubble is not well understood. Nevertheless, its pres¬ 
ence is found to affect pressure distribution over the sub¬ 
strate surface. The third region, known as the radial wall jet, 
includes the area outside the impingement zone which con¬ 
tains the jet flow, redirected laterally outward after impinge¬ 
ment. As the supersonic flow in the primary jet approaches 



Fig. 5. Schematic of a supersonic impinging jet flow field. 


the substrate, it decelerates through the formation of a bow 
shock. If the jet is not ideally expanded, oblique shocks in 
the jet plume (denoted as “jet shocks” in Fig. 5) interact 
with the bow shock resulting in the so-called “triple-shock 
structure”, with the third shock denoted as a “tail shock” in 
Fig. 5. The interaction between the three shocks is believed 
to govern the flowfield in the impingement zone [14]. 

A computational analysis of the impinging supersonic jet 
entails finding the solution to the following governing dif¬ 
ferential equations: 

• Reynolds averaged Navier-Stokes equations which in¬ 
clude the momentum conservation equations and a conti¬ 
nuity equation for compressible turbulent fluid low; 

• An energy conservation equation; 

• Closure relations which enable evaluation of the velo¬ 
city-fluctuations based terms (called the Reynolds stress 
tensor) in the Reynolds averaged Navier-Stokes and the 
energy conservation equations. Among various models for 
the closure relationship, one-equation Spalart model [15], 
and two-equation model (e.g. the k-e model [16], the k-co 
model [17], and the shear-stress transport (SST) model 
[18]) are most frequently used; and 

• Constitutive relations such as the ideal-gas law, the 
Sutherland’s law for the temperature-dependence of the 
(laminar) gas viscosity, and the models for the molecular 
and the turbulent Prandtl number. 

In our ongoing work, FEMFAB general-purpose finite- 
element computer program [19] is used to analyze the 
flowfield in supersonic axisymmetric impinging jets. The 
objective of this work is to establish relationships be¬ 
tween the primary-jet characteristics (more precisely the 
carrier-gas characteristics at the nozzle exit) and the flow- 
field in the stagnant zone behind the bow shock adjacent 
to the substrate. While a detail account of this work will 
be reported in a future correspondence, some of the results 
obtained will be utilized here to help construct an analytical 
function which relates the primary-jet and the geometrical 
parameters and the flowfield of a stagnant zone adjacent to 
the substrate surface. 

In accordance with many previous experimental and com¬ 
putational analyses, our work shows that the gas flow is sub¬ 
sonic in the stagnant region and that the component of the 
gas velocity normal to the substrate surface is quite small 
in comparison with the gas and particle exit-velocities. Un¬ 
der such condition Eq. (12), can be integrated to yield the 
particle impact velocity as [20]: 

4 mpact = u p e- 3 ^ LS, /%°r (25) 

where p st is the average gas density in the stagnant region 
while L sl is thickness of this region. In our ongoing work, 
the following relationships were found for p st and L st . 

p st = o e (— 1.04 + 2.27M e - 0.21 M e 2 ) (26) 

L st = R e (0.97 - 0.02Me) 1 < M e < 5 (27) 
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where M e is the gas Mach number at the nozzle exit. Thus, 
Eq. (23), in conjunction with Eqs. (22)—(24), can be used 
to compute the impact particle velocity for a given set of 
(stagnation) gas, process and feed powder parameter. In 
addition, as suggested by Dykhuizen [20]. Eq. (25) in con¬ 
junction with. Eqs. (22)—(24), can be used to determine the 
optimum particle size distribution. That is, if particles are 
too large (and heavy), they will not be accelerated enough 
with the nozzle, but their deceleration within the stagnant 
region will be modest. When the particles are too small (and 
light), on the other hand, they will acquire a high exit veloc¬ 
ity but may be greatly decelerated within the stagnant zone. 
Consequently, the impact velocities of particles which are 
either too large and too heave or too small and too light will 
be relatively low and the maximum impact velocities will 
be attained by particles of an optimum size (and weight). 
This can be seen in Fig. 6 in which the effect of the particle 
size, the gas type and the feed-powder particle size on the 
particle impact velocity is shown. The results displayed in 
this figure show that: (a) the particle size which is associated 
with the largest impact velocity scales with an inverse of the 
feed-powder density; (b) the maximum particle impact ve¬ 
locities for the two feed powders (copper and aluminum) are 
very close to each other for either choice of the carrier gas 
(helium or air); and (c) helium gives rise to a substantially 
higher particle impact velocities. All these findings are con¬ 
sistent with the results of Dykhuizen [20]. Also displayed in 
Fig. 6 are the experimental results for the impact velocity of 
copper particles with 1 p,m diameter accelerated in air and 
with 20 p,m diameter accelerated in helium [21]. Due to the 
observed relatively good agreement between the experimen¬ 
tal results and the corresponding model prediction displayed 
in Fig. 6, one might conclude that one of the key simplifi¬ 
cations used in the development of the model involving the 



Fig. 6. Effect of the particle size, the feed-powder material and the carrier 
gas on the particle impact velocity. 


neglect of variations of the gas density and the drag coeffi¬ 
cient along the diverging section of the nozzle was justified. 

5. Conclusions 

Based on the results obtained in the present work, the 
following main conclusions can be drawn. (1) Using a 
nonlinear regression analysis and a numerical solution 
for the one-dimensional isentropic gas flow in a cold-gas 
dynamic-spray nozzle, a relatively simple function is de¬ 
fined which relates the gas velocity at the nozzle exit with 
the nozzle expansion ratio and the carrier gas stagnation 
properties. (2) Analytical solutions for the particle velocity 
at the nozzle exit in the limits of very short and very long 
nozzle lengths can be combined to derive an analytical 
function which relates the particle velocity at the nozzle exit 
for intermediate nozzle lengths with various gas, process 
and feed-powder parameters. (3) To compute the velocity 
at which particles impact the substrate surface, deceleration 
of the particles in a stagnant subsonic region adjacent to the 
substrate surface must be considered. 
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